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In this work we study the Lepton Flavour Violating semileptonic r — ► ^70(980) 
decay within the context of SUSY-Seesaw Models, where the MSSM spectrum is ex- 
tended by three right handed neutrinos and their SUSY partners, and where the seesaw 
mechanism is used to generate the neutrino masses. We estimate its decay rate when it 
proceeds via the Higgs mediated channel r — ► fiH* — > /x/q(980), where H refers to the 
CP-even MSSM Higgs bosons hP and H°, and the Lepton Flavour Violating t^jlH ver- 
tex is radiatively generated via SUSY loops. In order to describe the /o(980) meson we 
follow the guidelines from chiral constraints. As an implication of our computation, we 
explore the sensitivity to the Higgs sector in this decay and compare it with other LFV 
tau decay channels. The confrontation of our predictions for BR(r — > /i/o(980)) with 
its very competitive present experimental bound leads us to extract some interesting 
restrictions on the most relevant model parameters, particularly, tan/3 and m H o. 
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1 Introduction 



The study of Lepton Flavour Violating (LFV) processes provides one of the most effi- 
cient indirect tests of supersymmetric (SUSY) extensions of the Standard Model of Particle 
Physics [1-6]. The reason is because in SUSY models the lepton and slepton mass matrices 
are not diagonal in flavour simultaneously, and this misalignment leads to intergenerational 
interactions between leptons and sleptons with neutralinos and charginos at tree level, that 
when placed into the loops of lepton flavour changing processes, can generate large rates. 
Furthermore, in the case of SUSY-Seesaw models, with extended lepton and slepton sectors 
by three right handed neutrinos, i/r, and their SUSY partners, Or, and where the seesaw 
mechanism is used to generate the neutrino masses (i.e., the so called Seesaw models of 
type I [7]), the size of the off-diagonal (in flavour) slepton mass matrix elements that are 
responsible for LFV, is governed by the strength of the neutrino Yukawa couplings which 
can be Y v ~ 0(1) or even larger for heavy M VR ~ 10 14 — 10 15 GeV. Thus, an interesting 
connection between neutrino and LFV physics follows, because the large Yukawa couplings 
of the Majorana neutrinos induce, via loops of SUSY particles, important contributions to 
LFV processes. In fact, these contributions are in some cases [1-6,8-23], at the reach of the 
present experimental sensitivity [24]. 

The LFV process that is the most sensitive to the neutrino Yukawa couplings, in the 
SUSY-Seesaw context, is fi — » cy, where the present experimental sensitivity is at 1.2 x 
10~ n [25, 26]. Also ix — e conversion in heavy nuclei, with present bounds at CR(/x — e, Ti) < 
4.3 x lCT 12 [27] and CR(/x - e, Au) < 7 x 1(T 13 ) [28], and /i -> 3e with BR(/i -> 3e) < 
1.0 x 10~ 12 [29], are quite sensitive to LFV in the ii — e sector. The most competitive LFV 
process in the r — ii sector is r — > /ry, whose upper bound is now set to 1.6 x 1(U 8 [30-32]. 
Moreover, the sensitivity to LFV in r — > 3/z has improved remarkably in the last years. The 
present upper bounds from BELLE and BABAR collaborations are 3.2 x 10~ 8 and 5.3 x 10~ 8 , 
respectively [33,34]. In the last years, several interesting bounds at the 10~ 8 level for some 
LFV semileptonic tau decays have also been provided [35-37]. 

In this work, we study the LFV semileptonic tau decay channel r — > /i/ (980) , which 
is competitive with other LFV tau decays due to the recently reported bound by BELLE 
collaboration [38], BR(r -> /// (980)) x £.R(/ (980) -> tt+tt-) < 3.4 x 10~ 8 . In fact, it 
is at present, the best bound in semileptonic LFV tau decays, improving the other present 
competitive bound of BR(t — > firj) < 5 x 10~ 8 [32]. The advantage of r — > fir] [9,13,21] 
and r — > /x/o(980) [18] over the r — > fi'-f channel is their potential sensitivity to the Higgs 
sector. Whereas the r — > /irj can be mediated by a Z boson and a CP-odd Higgs boson 
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A , and it is dominated by the A just at large tan (3 > 20 [21, 23], the r — > /// (980) decay 
is exclusively mediated by the exchange of the neutral CP-even Higgs bosons H° and h°. 
Therefore, through the r — > /x/ (980) channel one is testing directly the neutral CP-even 
Higgs sector at all tan (3 values . 

Our computation of the BR(r — > /i/ (980)) improves the estimate of [18] in several 
aspects. First, we demand compatibility with present data on light neutrino masses and 
mixings. Second, we do not use the mass insertion approximation, we take into account 
the full set of SUSY one-loop diagrams in the LFV vertex TfiH (H = h°,H°), and include 
the two contributions mediated by the ho and H respectively. Consequently, we explore 
the full 5 < tan/5 < 60 interval. Besides, the hadronization of quark bilinears into the 
/o(980) meson is performed here quite differently than in [18], where a simplified quark- 
flavour scheme was used to express these bilinears in terms of phenomenological meson decay 
constants. We instead pay close attention to the chiral constraints, following the standard 
Chiral Perturbation Theory (xPT) [39-41] and the Resonance Chiral Theory (R%T) [42-46] 
to incorporate resonances. Concretely, we follow the description of /o(980) in [45], where it 
is defined by a mixing between the octet and singlet components of the nonet of the scalar 
resonances which are included in R%T. Furthermore, we do not work in a generic Minimal 
Supersymmetric Standard Model (MSSM) but in constrained models with input parameters 
set at the high energies. Concretely we focus on two particular constrained SUSY scenarios 
of remarkable interest: the usual constrained MSSM (CMSSM) scenario [47], with universal 
soft SUSY masses at the gauge coupling unification scale, and the so-called Non-Universal 
Higgs Mass (NUHM) scenario [48] , with all the scalar soft masses being universal except for 
the Higgs sector ones. In this later case the physical Higgs boson masses, m^o and m#o, can 
be both light, ~ 100 — 250 GeV, indeed close to their present experimental lower bounds and, 
therefore, the corresponding Higgs mediated contribution to the previous LFV processes can 
be relevant, even for large soft SUSY masses at ~ 0(1 TeV). This is precisely the main 
interest of the channel r — > /i/o(980), namely, the fact that the decay rates can be sizeable 
even for large SUSY masses, Msusy ~ 0(1 TeV), in clear contrast with other competitive 
tau flavour violating channels like r — > /ry, whose rates decrease as 1/M| USY and lay below 
the present experimental bound for such a heavy SUSY spectrum. 



2 Framework for the r —> ///o(980) decay 

For the present study of the r — > /x/ (980) decay, we choose a SUSY-Seesaw framework where 
the spectrum of the MSSM is enlarged by three right-handed neutrinos, u Ri (i = 1, 2, 3), and 
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their SUSY partners, (i — 1,2,3). Here we assume a seesaw mechanism for neutrino mass 
generation and use, in particular, the parameterization proposed in [6] where the solution to 
the seesaw equation, relating the parameters of the six physical (mass eigenstates) Majorana 
neutrinos, Vi , and iVj (i — 1, 2, 3) to the neutrino Yukawa couplings, is written as 

m D =Y u v 2 = *y / ^i?V / ^ f/ PMNS- (1) 

Here, the Dirac mass, mo, the Yukawa neutrino coupling, Y v , and R are 3x3 matrices with 
full structure in flavour space. The orthogonal matrix R is defined by three complex angles 
9i (i = 1,2,3) [6]. m^ iag = diag {m vi , m V2 , m V3 ) denotes the three light neutrino masses, and 
^diag _ ^ mN ^ mA r 2) rriN 3 ) the three heavy ones. The two Higgs vacuum expectation 
values are t>i( 2 ) = v cos(sin)/3, with v = 174 GeV. The Pontecorvo-Maki-Nakagawa-Sakata 
unitary matrix £/pmns [49,50] is given by the three (light) neutrino mixing angles #12, #23 
and #i3, and three phases, 5, 0i and 02- With this parameterization is easy to accommodate 
the neutrino data. It further allows for large Yukawa couplings Y v ~ (9(1) by choosing large 
entries in m^ ag and/or Q; t . 

For the numerical predictions in this work we will set: 

m\ ~ , ml 2 = Ami, = 8 x 10~ 5 eV 2 , m 2 3 = Am 2 tm = 2.5 x 10~ 3 eV 2 , 

#12=30°, #23 = 45°, ^13 = 5°, 5 = 0! =02 = 0, (2) 

which are compatible with present neutrino data [24], and consider the two possibilities for 
the heavy neutrinos: 1) Degenerate, with m Nl = m N2 = m N3 = m N ; and 2) Hierarchical, 
with tojvi << TfiN 2 << ttin 3 - This later case is well known to provide a plausible scenario 
for the Baryon Asymmetry of the Universe (BAU) via leptogenesis. 

Regarding the SUSY parameters we will work within two different constrained MSSM- 
Seesaw scenarios, the CMSSM with universal soft SUSY breaking parameters (including 
the extended sneutrino sector) and the NUHM model with non-universal Higgs soft masses. 
Thus, in addition to the previous neutrino parameters, m N . and 0j, the input parameters of 
these two models are respectively, 

CMSSM : M , M 1/2 , A , tan (3 , sign(^) , 
NUHM : M , M 1/2 , A , tan/3 , sign(/i) , M 2 Hi = M 2 (1 + 5J, M 2 H2 = M 2 (1 + 5 2 ). (3) 

where M , Mi/ 2 and A are the universal soft SUSY breaking scalar masses, gaugino masses 
and trilinear couplings, respectively, at the gauge coupling unification scale, M x — 2 x 10 16 
GeV. The other parameters are, as usual, the ratio of the two Higgs vacuum expectation 
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values, tan/5 = v 2 /vi, and the sign of the \i parameter, sign(/i). Notice, that the departure 
from universality in the soft Higgs masses of the NUHM is parameterized here in terms of the 
two dimensionless parameters 5i and 62 ■ Consequently, by taking 5\ = 62 = in the NUHM 
one recovers the CMSSM case. Finally, in order to evaluate the previous SUSY parameters 
and the physical masses at low energies (taken here as the Z gauge boson mass mz), we solve 
the full one-loop Renormalization Group Equations (RGEs) including the extended neutrino 
and sneutrino sectors. For this and the computation of the full spectra at the low energy 
we use here the public FORTRAN code SPheno [51]. In the numerical estimates we will set 
M = Mi/2, Aq = and sign(/i) = +1, for simplicity. 

For the purpose of the present analysis the most relevant difference between the two pre- 
vious constrained SUSY-Seesaw scenarios is the spectrum of the Higgs sector. In particular, 
we want to explore the interesting case where the neutral Higgs bosons that mediate the 
r — > /i/ (980) decay are light, while keeping the SUSY spectra heavy enough as to suppress 
the other competitive LFV tau decay channels like, for instance, r — > fij. This is clearly 
possible within the NUHM-Seesaw scenario, as illustrated in Fig. 1. We see in this figure 
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Figure 1: CP-even Higgs boson masses in the NUHM-Seesaw scenario: 1) mjyo as a function of Mstjsy = 
Mo = M 1 / 2 for several input 5^2 (left panel). The predictions in the CMSSM-Seesaw scenario (Si =$2 = 0) 
are included for comparison; 2) m H o and m h a as functions of tan (3 for Msusy = 250 GeV and 750 GeV 
(right panel). 



that, by properly adjusting the input 5\ and 62 parameters, the heavy Higgs boson H° can 
get masses as low as 100-250 GeV even for a very heavy SUSY spectrum. For instance, 
for 6 1 = -2.4, 5 2 = 0, tan/5 = 50, M SUSY = M = M 1/2 = 750 GeV and the other input 
parameter values as specified in this figure, we get m#o = 249 GeV and m^o = 122 GeV, to 
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be compared with m H o = 998 GeV and m h a = 122 GeV of the CMSSM-Seesaw case. With 
other specific choices for 5 2 ^ one gets even lower values of m H o [21]. For the following 
numerical analysis and, for simplicity, we will set, however, S 2 = and play just with Si. It is 
worth also mentioning that the predictions for m A o (not shown in this figure) are practically 
indistinguishable from those of m H o [21]. 

Within the previous scenarios for the neutrino and SUSY sectors, it is well known that 
one can get large LFV decay rates if one chooses large entries in m% ag and/or complex 9i, 
basically due to the large size of Y v in these models. This can be understood more easily 
in the Leading Logarithmic (LLog) approximation where, the tau-muon flavour violation, 
which is of our interest here, is qualitatively well described by the parameter, 

fa = - ^ M l +Al) (*?^)32, ^ log (^WM = 1,2,3, (4) 

where Msusy is an average SUSY mass. The size of | ^32 1 can be indeed quite large. For 
instance, for mass values of the heavy neutrinos rriN 3 (or mjv) in the range 10 14 — 10 15 GeV 
and 9i (i= 1 or/and 2) with large modulus in the range 3 — 5 or/and large argument in the 
range [±7r/4, ±7r/2] one can get values of \S 32 \ as large as 0.5-10. This is clearly illustrated in 
the contour plots of Fig. 2, where we have considered both scenarios with either degenerate 
or hierarchical heavy neutrinos and we have explored in the (m^,^) parameter space. In 
the hierarchical case the relevant mass is the heaviest one m^ 3 and the predictions for \8 32 \ 
do not vary appreciably with mjv lj2 - I* 1 addition, we have checked that \S 32 \ is nearly constant 
with 9 3 . The contour plots for Q x (not shown) are very similar to those of 9 2 . We have also 
found that the largest values of \S 32 \ are obtained for the degenerate case with both 61 and 9 2 
being large. This is also clearly illustrated in the lower right panel of Fig. 2. For instance, we 
get \S 32 \ c± 5 for m N = 10 14 GeV and 9 X = 9 2 = 3exp (in/A). Notice also that values of \S 32 \ 
larger than ~ 0.5 correspond in our parameterization of the Yukawa coupling matrices in (1) 
to values of |Y"^| 2 /(47t) that are above the threshold where the SPheno code sets the limit of 
perturbativity, which is at \Y v \ 2 /(Ait) ~ 1.5. It means that, in the following, we will be able 
to provide full predictions for the decay rates with the SPheno code only for those model 
parameters producing Y v values that are within the perturbativity region or, equivalently, 
leading to \S 32 \ < 0.5. The implications for the r — > /x/ (980) decay of values |<5 32 | > 0.5 will 
be explored later, not with our full computation implemented by us in SPheno, but using an 
approximate formula that will also be presented here and that turns out to work reasonably 
well. 

Next, we specify our framework for the hadronization of the quark bilinears into the 
/o(980) meson. We use here the chiral Lagrangian of R%T that is a suitable tool to realise 
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Figure 2: Contours of | J32 1 in the CMSSM-Seesaw scenario: 1) For hierarchical heavy neutrinos. Upper left 
panel: in the (|#2|>TOjv 3 ) plane for arg# 2 = tt/4. Lower left panel: in the (arg0 2 : m Af 3 ) plane for \9 2 \ = 3. 
The other heavy neutrino parameters are set to 8\ = 83 = 0, mjVi = 10 10 GeV, m,M 2 = 10 11 GeV; 2) For 
degenerate heavy neutrinos. Upper right panel: in the (|^2|,^w) plane for arg#2 = 7r/4 and 6\ = 83 = 0. 
Lower right panel: in the (|#i| = \82\,itin) plane for arg#i = arg6*2 = 7r/4, and ^3 = 0. In all plots we have 
set: Msusy = Mq = Mx/%, Aq = 0, tan/3 = 50, and the 9% are expressed in radians. 



7 



the l/N c expansion of SU(Nc) QCD and includes both the Goldstone bosons $(7r, K and 
i]) and the resonances as active degrees of freedom, and their interactions. For the present 
work, it is sufficient to consider the lightest nonet of scalar resonances R(0 + ) in R%T, 



-R X T - + ^kin + £(2) > ( 5 ) 



where, 



F 2 

4 2) = "J" + X+) , F~F„~ 92.4 MeV, 

^ki„ = liWRV.R-MlR 2 ) , 

£f 2) = c d (i2«X) + c m (i?x+) , (6) 
and (...) is short for a trace in the flavour space. The other quantities in (6) are: 



i[u\d^ - ir^)u - w(<9 M - ^)w f ] , u = exp [i $/ (a/2F) 

X+ = w'W + u^it , X = 25 (s + ip) , 



V M it> = d^R+lT^R], r M = ^[«t(^-ir> + «(^-^>t] > ( 7 ) 

being $ the nonet of Goldstone bosons, £ M = i> M — a M , r M = t> ^ + a M , and t> M , a M , s and p 
are the nonets of vector, axial-vector, scalar and pseudoscalar external fields, respectively. 
Short-distance dynamics [44] constraints the couplings of R%T by imposing the QCD ruled 
behaviour of Green functions and associated form factors. For the couplings in £^ one 
gets 1 : 

2c m = 2c d = F . (8) 

Finally the chiral tensor \ gives masses to the Goldstone bosons through the external scalar 
field. In the isospin limit one has : 

2B m u = 2 B Q m d = m 2 n , 

2 B m s = 2m 2 K - ml. (9) 

The QCD spectrum of scalar resonance states is far from being settled and constitutes, at 
present, a highly debated issue. It is not our goal in this article to enter in the details of the 
discussion and, therefore, we will attach to the scheme put forward in [45] for the description 



1 Short-distance constraints on the RxT couplings depend on the operators included. The result in (8) 
is obtained when only linear operators in the resonances are considered [42]. A weaker constraint, though 
compatible with that result, arises if non-linear couplings in the resonances are included [46]. 
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of the isosinglet / (980) state. The later is defined as a rotation of the octet R$ and the 
singlet R components of the -R(0 + ) nonet, 




cos 6 S sin 6 s 
— sin 6 s cos 6 s 



/o(1500) 
/o(980) 



(10) 



The value of the 6s mixing angle is uncertain. In the analysis carried out in [45] considering 
nonet breaking (i.e. subleading effects in the large-A^ expansion) a possible dual scenario 
is favoured : 



A) The candidates for the nonet are: / (980), X o *(1430), a (1450) and / (1500). In this 
framework the a (980) is dynamically generated (through loops). The mixing angle, 
around 6 S ^ 30°, provides a dominant non-strange component for the / (980) state 
and, consequently, justifies its dominant decay into two pions. 

B) The nonet would be composed by: / (980), a (980), ifj(1430) and / (1500). Hence 
do (980) is a pre-existing state in the Nq — > oo limit. The mixing angle in this case is 
around 6s — 7°, that gives a noticeable strange component for the /o(980) state. 

Given the uncertainty provided by the large corrections due to l/N c subleading effects we will 
consider the two previous scenarios for the /o(980) as plausible and will present estimates 
of the r — > /// (980) decay rates for the two mixing angles, 6 S ^ 7° and 6 S ^ 30°. The 
dispersion between these two results can be considered as part of the theoretical error in our 
estimates. 

Finally, the hadronization of the relevant scalar quark bilinears into the / (980) is imple- 
mented by replacing the following expressions in the results for the decay rates at the quark 
level, 



uu 



dd 



s s 



with 



Si 



T2 BoC 



m Ri i 



0,3,8, 



(12) 



and, according to (8), c m = F/2. As i? 3 does not contain information on / (980) (in the 
isospin limit) we will discard the S3 contribution. 
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Before proceeding a word of caution is necessary when dealing with processes with reso- 
nances as initial or final states. A resonance is not an asymptotic state as it decays strongly. 
Hence from a quantum field theory point of view R/\T only describes the creation, propaga- 
tion and destruction of resonances and the later should not appear as "in" or "out" states. 
For instance, in our case the physical process should be r — > finn mediated by a / (980) 
state, and not r — > /x/ (980). Then it would proceed to study the scalar state as was done 
with the vector ones in [21]. However the description of scalars, as has been pointed out, 
is far from clear and therefore considering the /o(980) as an asymptotic state should not 
increase effectively the already rather large uncertainty. 



Analytical results 

The semileptonic r — > /i/o(980) decay can be mediated by h° and H° Higgs bosons, as 
shown in Fig. 3. In this figure the LFV vertex is represented by a black circle and the 
hadronic vertex by a grey box. The total amplitude for this decay, T# = T h o + T H o, is 
first evaluated at the quark level, that is for r — > fjiqq, and then at the hadron level by 
substituting the quark bilinears by the corresponding scalar currents containing the /o(980) 
meson as evaluated from Cr x t in (5). The amplitude at the quark level can be computed 
in terms of the corresponding t/j,H p one-loop vertex functions, H^ R , with H p = h°, H°, 
resulting from the evaluation of the diagrams in Fig. 4 with sleptons, lx, sneutrinos, Ux, 
charginos, Xai an d neutralinos, Xa, m the loops. The resulting amplitude at the quark level 
is given by: 



Th(t - m) = E K W S£ Wlt] [qP L q] + H^S® [HP R T] [qP R q] 

III TT ^ 

h°,H° h p 

+ HPS® [PP L t] [qP R q] + lf«S« [pP R r] [qP L q] } . (13) 



3 Results for BR(r 



/x/o(980)) 



where P L)R = (1 =F 7s)/2, and 




(14) 
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Figure 3: Higgs-mediated contributions to the LFV semileptonic r — > ^/ (980) decay 




Figure 4: Relevant SUSY one-loop diagrams for the Higgs-mediated contributions to the r — > /x/o(980) 
decay. Here H p = h°,H°. 
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sin a 




a {p) _ 
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cos a > 




— cos a 






sin a 


V 


i sin /3 


J 




V 


—i cos (3 J 



Here mw is the W gauge boson mass, m q is the q quark mass, a is the mixing angle in the 
Higgs sector, and g is the SU(2) gauge coupling. The three entries in are ; i n order from 
top to bottom, for H p = h°, H°, A , respectively. 

The results of the LFV vertex functions are taken from [14], and are not written here 
explicitely for shortness. Just to mention that it is a full one-loop computation, including all 
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the contributions with charginos in the loops, Hl{r) c > anc ^ those with neutralinos, Hl(r) n- 
Besides, all these contributions are written in terms of the physical particle masses. As 
we have mentioned before, these physical masses are computed here in the SUSY-seesaw 
scenario by solving the one-loop RGEs with SPheno and for a given set of universal (in the 
CMSSM) or non-universal conditions (for the NUHM) at the unification scale. Since the 
three right-handed neutrinos and their SUSY partners are included in the RGEs, they will 
affect as well in the predicted physical masses at the low energies. 

To get the amplitude for the process r — > /x/ (980) we substitute the quark bilinears of 
(11) in (13) and use (10) and (12). Notice that it is just the scalar part in \^Pl,rq], and not 
the pseudoscalar, the one that contributes in the present case. We obtain: 



T h {t - /x/ (980)) 



p=h°,H° 



where 



and 
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2mw 2m 2 H 



(^ + ^)( 



H® + HP 



)■ 



(16) 



(17) 



J 



(p) 



a. 



(p)* 



y/3 sin/3 



V2 



sin 9 S + cos 9 S 



a 



(p)* 



cos {3 



3 

7! 



sin 9s m 2 + ^cos^s — v^sin^sj 2 



in 



K 



t(p) _ 7 

J R — J 



(p)* 



(18) 



Notice that due to the mass relations in (9), the couplings of the Higgs bosons, h° and H Q , 



to the quarks (q — u,d,s ), S^ q and S K Lq 



; / ' 1 in (14), being proportional to the quark masses, 
lead to Higgs- /o couplings that are proportional to m 2 P (P = ir,K). This is seen clearly 
in the predicted functions J^ R of (18). In consequence, the dominant contributions to 
BR(r — > /i/ (980)) will come clearly from the terms in the amplitude that are proportional 
to m\. 

Finally, the result of the branching ratio for the r — > /i/ (980) decay is given by, 



1 X 1/2 (m 2 m 2 m 2 f ) 1 ^ 
BR(^,,/o(980)) = - 4^ fJ ~ 2 j:\T H \\ 



m 2 T T T 



where 



(to m + m r ) 2 - m 



i,/ 



4m T 



■fo u i „ |2 
— |C/jO + C H o\ 



(19) 



(20) 
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T r is the total r width and X(x, y, z) = (x + y — z) 2 — 4xy. 

Approximate formula 

Next we derive a simple formula which approximates reasonably well our full one-loop 
prediction in (19) and (20). For this, we work within the approximation of large tan/3 that 
is appropriate for LFV tau decays, whose rates grow quite fast with this parameter. This is 
especially relevant for channels where the LFV rates are dominated by the Higgs mediated 
diagrams, as it is the present case, and where the growth with tan f3 is extremely pronounced. 

The other approximation which is used frequently in the literature, due to its simplicity, 
is the use of the mass insertion (MI) method, where the tau-muon LFV is encoded in the 
dimesionless parameters 5^ 2 Y (XY = LL, RR, LR). In the SUSY models the dominant one 
is 5^2 an d its expression in the LLog approximation, (5f 2 L )LLo g = 832, is that given in (4). 

It is known [14] [16] that at large tan f3 the vertex function dominates Hr by about 
a factor m r /m M . In addition is by far larger than in this limit, and one can safely 
neglect the later one. More specifically, by using the MI approximation, its chargino and 
neutralino contributions in the large tan/5 and heavy M SU sy limits give, correspondingly, 
the following expressions : 

HiV = ^^5 32 tan 2 /?, 
L ' c 16tt 2 12m w 32 H ' 

<' = I&2^( 1 - 3t - 2 "»')^ t - 2 ' 5 - < 21 > 
One can further verify that H c dominates H n by about a factor 20, so that we will simplify 
H L ~ H LtC . 

On the other hand, we also consider the large tan (3 limit of the functions that define the 
H° couplings to / (980), Jl and Jr in (18). We obtain : 



J 



= J 



(HO) 



R 



2v/3 



tan/? 



y/2 



sin 9 s ml + (cos 6 S — V2sm6s)2m 2 K 



(22) 



By using the above sequence of approximations and by neglecting the muon mass, we 
finally get the following simple result: 



BR(r -> /i/ (980)) 



approx 



1 / 2 2 \2 

16wK- m /o) 



9 



2m\y m H o 



L_ j(H°) tt(H°) 
2 J L n Lfi 



1 



(23) 



7.3 x 10- 8 (6 S = 7°) 
4.2 x 10~ 9 (9 S = 30°) 



J32 1 



100 



m^o(GeV) J V 60 



4 / tan/3\ 6 
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In the last line we see explicitly the fast growth with tan/3, as (tan/3) 6 , the expected de- 
pendence with the relevant Higgs mass, as (m^o) -4 , and also with the LFV parameter, as 
| <532 1 2 - The two numerical factors correspond to the two assumed values for the mixing angle 
that defines the / (980) state, 9 S = 7° and 9 S = 30°. These two results differ by a factor 17, 
meaning that the predicted rates will carry a theoretical uncertainty of about this number, 
due to the uncertainty in the definition of the / (980) state. 

Numerical results 

In the following we present the numerical predictions for BR(r — > /i/ (980)). We first 
show the results from the full computation in (19) and (20) and then compare with the 
approximate results in (23) and also with the rates of other LFV tau decay channels. 

In Fig. 5 it is shown the BR(r — * /i/o(980)) versus the heavy neutrino masses, in both 
scenarios with hierarchical and degenerate heavy neutrinos. In the hierarchical case we 
display just the dependence with the relevant mass, m^- As expected, from the previously 
manifested behaviour of | (5 32 1 with m^. A (or with m^, in the degenerate case) in Fig. 2, we 
find a fast growing of BR(r — > ///o(980)) with this mass. Although not explicitely shown 
here, we have also checked in the hierarchical case, the near independence on the other 
masses, m Nl and m N , 2 . From this figure it is also evident that by choosing properly the 5i 
and S 2 parameters of the NUHM scenario, such that the relevant Higgs boson mass m H o gets 
lower than for S± — 8 2 — 0, the branching ratios get larger than in the CMSSM scenario. 
Finally, by comparing the rates of the two neutrino scenarios, and for the same input model 
parameter values, including the same m N and m N3 , we find rates in the degenerate case 
that are generally larger than in the hierarchical case. For instance, for the choice of input 
parameters in Fig. 5 we find larger rates by a factor of about 3. In the following we will 
focus more on the hierarchical case since it has the appealing feature of providing successful 
baryogenesis, via leptogenesis, for some regions of the heavy neutrinos parameter space. 

We present the predictions of the BR(r — > /i/ (980)) versus de soft SUSY masses M 
and M1/2 m Fig. 6. Here we take again M = My 2 = -Wsusy and compare the results in 
both scenarios, the NUHM with S 1 = —2.4 and 5 2 = 0, where the predicted Higgs boson 
masses for large tan/3 ~ 50 lay within the interval 100-250 GeV, and the CMSSM. The 
most evident feature in this plot is the different behaviour of the BR(r — > /i/ (980)) with 
Msusy in these two scenarios. Whereas in the CMSSM the rates are found to decrease with 
increasing M SUSY , as expected, it clearly does not happen in the NUHM. In fact, the rates 
are practically constant for M SUSY > 400 GeV. The reason for this behaviour is that the 
SUSY particles do not decouple at large M SUSY in this decay. The non-decoupling behaviour 
can be checked analytically in that the LFV vertex, described by the dominant form factor 
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Figure 5: BR(r — > /x/ (980)) in the NUHM-Seesaw, for several Si values, and in the CMSSM-Seesaw versus 
the relevant heavy neutrino mass, 1) for hierarchical heavy neutrinos (left panel), and 2) degenerate heavy 
neutrinos (right panel). 

Hl, tends to a constant value at asymptotically large Msusy, as indicated in (21). Since, 
on the other hand, m#o is kept at the low region even for large Msusy, then a constant Hl 
with Msusy implies approximately constant BR(r — > /x/ (980)) as well. 

Another interesting feature of the predicted rates in the NUHM scenario, that is man- 
ifested in Fig. 6 as well, is the clear dominance by many orders of magnitude of the H° 
contribution over the h° one in the whole Msusy considered interval. This is due to the 
fact that at large tan (3 the H° contribution is enhanced by a tan 6 (3 factor, whereas the 
h° one is suppressed in this limit. In fact, we also see in this plot that the total rates are 
nearly indistinguishable from the H° contributions. Thus, to neglect the h° contribution is 
an extremely good approximation. 

Concerning the Higgs sector parameters, the BR(r — ► /i/o(980)) is mainly sensitive to 
tan/3 and m#o since, as said before, the if°-mediated LFV semileptonic decays grow very 
fast with both tan/? and l/m#o. In fact, in the approximation given in (23), as already said, 
BR(r — > yu/ (980)) goes as (tan/?) 6 and (1/m^o) 4 , respectively. 

The predictions of BR(r — ► /i/ (980)) as a function of tan/5 are shown in the right 
panel of Fig.6. We show again separately the h° and H° contributions and the total rates 
which are clearly dominated by the H° in the full studied interval of tan (3. Besides, it also 
displays the fast growing of the total rates with tan/3 , reaching values at the ~ 10~ 9 level 
for tan (3 ~ 50 which are close but still below the present experimental bound. We also see 



15 



10' 



1(T 



10 



10 



•^io- 11 

t 10" 12 

co 

10- 14 
io- 15 



10" 



10"' 



8, =-2.4, Full a 

8, =-2.4, H° • 

8, =-2.4, h° + 

CMSSM o 



m N = (10 10 , 10", 10"*) GeV 
A = 0, tan p = 50 
e 2 =3e i " /4 , e 1 =e 3 =o 
5 Z = o, e s =7° 



,1 1 



J °o , 



°oo 



°°ooo ( 



°°ooo, 



'°ooo Q 



°oooooi 



++- 



ooooooooo 



++++ ++-H^ H+H+ ^ |+|+ ^^ 



300 



400 



500 



600 



700 



800 



900 



M = M 1/2 (GeV) 



10"' 



10"' 



3>W 
I10" 14 

00 



io- 



10- 



■ m N = 


(10 10 , 


10 11 , 10 u ) GeV 




M = 


M 1/2 = 


250 GeV, A = 






3e , 


6, = e 3 = _ 




\ 


-2.4, S 2 


= o, e 3 s = 7° 



























/ 








Full • 


* 






Approx. a 


* 






h° 


* 






H° x 




10 


20 30 


40 50 






tan B 





Figure 6: BR(r ->■ /x/ (980)) in the NUHM-Seesaw scenario: 1) As a function of M = M 1/2 = Af SU SY (left 
panel). We show separately the Hq and ho contributions as well as the total. The predictions for the total 
rates within the CMSSM-Seesaw scenario are also included for comparison; 2) As a function of tan (right 
panel). Again, the dominant H , the subdominant h° and the total rates are displayed. We also include 
here the approximate predictions given by (23) for comparison with the full rates. The dotted horizontal 
line at the top of the graphics is the present experimental upper bound. 



that the particular shape of the curve for the total rates is a consequence as well of the m#o 
dependence with tan/? in these SUSY scenarios, as illustrated in Fig. 1. 

The comparison between our predictions for the full result in (19) and (20) and the 
approximate result in (23), which includes just the H° boson contribution, can be seen 
as well in Fig. 6. The agreement between the full and the approximate results is quite 
remarkable, for all the studied values in the 5 < tan (3 < 50 range. Therefore, we conclude 
that our simple formula (23) provides a very good approximation to BR(r — ► /i/ (980)) for 
all tan f3. 

It is interesting to compare r — > yu/ (980) to other Higgs-mediated LFV tau decay chan- 
nels like r — > fit] and t —>■ 3fi. First, notice that our previous result of the H° dominance in 
the r — > /x/o(980) channel over the full tan/? interval, is not true for the correlated channel 
r — > fit], nor the leptonic r — ► 3/i decay. The semileptonic LFV r — > fir] decay can be 
mediated by a CP-odd A Higgs boson and a Z boson, but the contribution from A dom- 
inates the full rates only in the large tan/? > 20 region [21,23]. The r —>■ 3fi channel can 
be mediated (apart from the box diagrams, which are negligible) by a photon, a Z boson 
and the three neutral Higgs bosons, h°, H° and A [16]. The photon dominates largely this 
decay, except at the extreme high values of tan j3 > 60 and Msusy > 1 TeV, where the two 
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type of contributions from the photon and the Higgs bosons, H° and A compete. These 
features can be seen clearly by comparing the corresponding approximate formulas, valid at 
large tan /5, for their respective Higgs boson contributions. That is, one should compare our 
result in (23) to the previous results of BR(r -> firj) [13,21] and BR(r -> 3/x) [8, 10, 12, 16] 
for the same input parameters. These are [21], 



BR(r -> nrj) Hti 
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From this comparison, we conclude that, for the same choice of the model parameters, and for 
9 S = 7°, the three rates BR(r — > /i/ (980)), BR(r — > ^77) and BR(r — > 3/i) are very similar 
if tan/5 > 60 and M 5C/5Y > 1 TeV. Concretely, we predict BR(r -> /// (980)):BR(t -> 
3/x):BR(r — > ^77) ~ 0.6 : 1 : 1, and they are all at the ~ O(10~ 7 ) level for |<5 32 | ~ 1, 
m H ~ 100 GeV and tan/3 ~ 60. Therefore, the three are closely competitive channels. It 
should also be mentioned that our estimate of BR(r — > /x/ (980)) for #5 ~ 7° and for the 
same input parameters, mn, tan/5 and jfe], is about one order of magnitud smaller than 
the prediction in [18]. They also predict a different ratio among the three LFV channels of 
~ 1.3 : 0.5 : 1. We believe that the main differences come from our different approaches for 
hadronization which produce, as we have already said, a dispersion in the results by a factor 
of 0(10). 

Finally, we summarize the sensitivity to the Higgs sector in the NUHM-Seesaw scenario 
in Fig. 7. In this plot we are using the approximate formula in (23) and we are setting 
#2 = 3 and Si = —2.4, 5 2 = 0. The soft masses are varied in the range 200 GeV < M = 
M1/2 = A^susy < 750 GeV. The explored m H o values in this plot correspond precisely to 
the output Higgs masses for this later M SU sy interval. The main conclusion from this plot is 
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Figure 7: Sensitivity to the Higgs Sector in r — > [i fo (980) within the NUHM-Seesaw scenario. The predicted 
rates are within the approximation of (23) and are displayed as a function of to#o, for various choices of 
large rriN 3 and tan/3. The dotted horizontal line is the present experimental upper bound. 



that for large m^ 3 ~5x 10 14 — 10 15 GeV and large tan (3 ~ 50 — 60 the predicted rates are 
already at the present experimental reach and, therefore, there is indeed Higgs sensitivity 
in this channel. In this concern, we find interesting to further explore if with the present 
experimental bound of BR(r -> j// (980))x BR(/ (980) -> < 3.4 x 10" 8 one may 

already exclude some region of the model parameter space. Our conclusion is that indeed 
it is possible to exclude the regions in the (m H o, tan/5) plane as summarized in Fig. 8. In 
this plot we assume, for simplicity, BR(/ (980) — ► -n + it~) ~ 1 and choose the specific input 
values, |5 32 | = 0.1,0.5,1,5,10. For each fixed |5 32 | the excluded region is the area above 
the corresponding contour line. For completeness, we have also included in this plot the 
present experimental lower bound for the SM Higgs mass at 114.4 GeV. Some words of 
caution should be said, anyway, about the conclusions from this plot since there are large 
uncertainties involved in the theoretical estimate of BR(r — > /i/o(980)). There are two main 
ones: 1) the uncertainty in the definition of /n(980) that, as evaluated in (23), can produce 
a dispersion of more than one order of magnitude in the predicted rates, and 2) the use of 
the approximate formula for values of jfe] > 0.5 which are out of the region that is allowed 
by a perturbative approach. The use of the MI approximation for such large values of 
is also questionable. 
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Figure 8: The excluded regions in the (mjjo, tan/3) plane are the areas above the contour lines corresponding 
to fixed |<5 32 | =0.1,0.5,1,5,10. 



4 Conclusions 

In this work we have studied in full detail the LFV semileptonic tau decay channel r — > 
/i/o(980) within the context of two constrained SUSY-Seesaw models, the CMSSM-Seesaw 
and the NUHM-Seesaw which have very different Higgs sector spectra. Concretely, we have 
selected NUHM-Seesaw scenarios with a light Higgs sector, h°, H° and A , in the 100-250 
GeV range, and considered several possibilities for the SUSY sector, varying the SUSY mass 
^susy in the 200-1000 GeV range. Through all this analysis, we have required compatibility 
with both the present experimental upper bound for this decay and with neutrino data for 
masses and oscillations. 

We have presented a full computation of BR(r — > /x/ (980)) that includes the complete 
one-loop SUSY diagrams with charginos, neutralinos, sleptons and sneutrinos contributing 
in the loops to the relevant effective LFV r/zif vertex. We have also taken into account the 
two kind of Higgs-mediated diagrams, with h° and H° in the internal propagator connecting 
the LFV vertex with the final quark-antiquark pairs. On the other hand, and in order to 



19 



provide predictions for the final meson / (980), we have performed the hadronization of the 
quark bilinears by means of the standard techniques in yPT and RyT. We have shown that 
in this chiral approach, the Higgs coupling to the /o(980) is dominated by its strange quark 
component. The leading term in this coupling is proportional to m 2 K) which is a consequence 
of the Gell-Mann-Oakes-Renner mass relation (B m s = m 2 K — 1 /2ml), and the fact that the 
Higgs coupling to the strange quark is proportional to m s . On the other hand, the H° — f 
coupling is dominant over the h° — fo coupling since the first one goes approximately as 
tan/5 in the large tan/3 limit (due again to this behaviour of the Hss coupling), whereas the 
second one is suppressed in this limit. 

We have analysed in full detail the dependence of BR(r — > /i/ (980)) with all the pa- 
rameters defining the two constrained SUSY-Seesaw scenarios and we have extracted from 
this analysis which are the relevant ones. Regarding the heavy neutrino sector, and for the 
most BAU favorable scenario of hierarchical heavy neutrinos, the most relevant parameters 
are the heaviest neutrino mass, m N3 and the 8 li2 angles. Concerning the SUSY and Higgs 
sectors the most relevant parameters are the SUSY masses, driven by Msusy, the CP-even 
Higgs boson mass, m H o, and tan/?. 

In the numerical predictions, we have found much larger rates in the NUHM-Seesaw 
than in the CMSSM-Seesaw scenarios, due mainly to the lighter Higgs mass m H o found 
in the first scheme. Indeed, it is just in the NUHM-Seesaw case where the predictions for 
BR(r — > /i/ (980)) can reach the present experimental sensitivity. We have shown, that in 
order to get values of BR(r — > yu/ (980)) at the 1CT 8 — 1CT 7 level one needs large values for the 
relevant parameters, namely, tojv 3 ~ 10 14 — 10 15 GeV, |#i j2 | ~ 2 — 3, ± arg(6*i j2 ) ~ tt/A — 3tt/A, 
tan /3 ~ 50 - 60 and m H o ~ 100 - 200 GeV. 

In addition to the full results, we have provided an approximate simple formula for 
BR(r — > /i/ (980)) which has been obtained in the large M S usy and large tan/3 limit, and 
with the MI approximation for the relevant LFV parameter <5 32 . Furthermore, we have shown 
in this work that this approximate result agrees pretty well with the full result in practically 
all the explored parameter space. The main basic features of the full predicted rates are 
very well reproduced by the simple formula in (23), which summarizes the fast growing 
with tan/3, going as (tan/3) 6 , with l/m H o, going as (1/m^o) 4 , and being approximately 
constant with M SU sy- The dependences with m N , A and # lj2 g° via the 5 32 parameter, and 
the large mjv 3 values are what enhance dominantly the rates, growing approximately as BR 
~ |mjv 3 logmArJ 2 . 

The most important conclusion from this work, as illustrated in Figs. 7 and 8, is that the 
LFV tau decay r — > /i/ (980) is indeed sensitive to the Higgs sector of the NUHM-Seesaw 
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models. Concretely, it is mostly sensitive to the CP-even Higgs boson H°, and therefore it 
complements very nicely the previous searches via the r — > fj,rj decay which is sensitive to 
the CP-odd Higgs boson A . These two channels together with the leptonic r — > 3/x decay 
are undoubtly the most competitive LFV tau decays where to look for indirect Higgs signals. 
As a final product of our analysis we have extracted some excluded areas in the parameter 
space of these models by using our approximate formula. The sensitivity found here to the 
Higgs sector will presumably improve in the future if the experimental reach increases up to 
1CT 9 — 1CT 10 , as it seems to be the case in the future SuperB and flavour factories [52]. 
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